A simple Monte Carlo model for crowd dynamics 
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In this paper we introduce a simple Monte Carlo method for simulating the dynamics of a crowd. 
Within our model a collection of hard-disk agents is subjected to a series of two-stage steps, implying 
(i) the displacement of one specific agent followed by (ii) a rearrangement of the rest of the group 
through a Monte Carlo dynamics. The rules for the combined steps are determined by the specific 
setting of the granular flow, so that our scheme should be easily adapted to describe crowd dynamics 
issues of many sorts, from stampedes in panic scenarios to organized flow around obstacles or through 
bottlenecks. We validate our scheme by computing the serving times statistics of a group of agents 
crowding to be served around a desk. In the case of a size homogeneous crowd, we recover intuitive 
results prompted by physical sense. However, as a further illustration of our theoretical framework, 
we show that heterogeneous systems display a less obvious behavior, as smaller agents feature shorter 
serving times. Finally, we analyze our results in the light of known properties of non-equilibrium 
hard-disk fluids and discuss general implications of our model. 



PACS numbers: 89.90.+n,02.70.Uu 

Keywords: Simulations of granular matter, crowd dynamics, 
disks 



social behavior, Monte Carlo methods, hard 



I. INTRODUCTION 



Crowd dynamics is the object of a comprehensive body 
of literature, from classic broad treatises covering gen- 
eral topics [1-4] to more specific applications, ranging 
from psycho-dynamical models [5] and control issues [6, 7] 
to panic scenarios [8, 9]. Likewise, several in-depth re- 
views discuss computational models of pedestrian dy- 
namics [10-13], ranging from lattice-gas models [10] to 
molecular-dynamics-based simulation schemes such as 
the behavioral (or social) force model [14, 15] and cel- 
lular automata methods [16]. 

It is the purpose of the present paper to introduce a 
novel model of crowd dynamics, based on a simple Monte 
Carlo (MC) scheme of hard disks dynamics. Our aim is 
two-fold: one the one side, we wish to outline the essen- 
tials of an original modeling strategy, that has never been 
considered so far to our knowledge. The simple philos- 
ophy behind our model can be effortlessly extended and 
adapted to many different scenarios, providing a viable 
and fast alternative to many current schemes for explor- 
ing real-life contexts of many sorts in crowd dynamics 
issues. At the same time, we also apply our theoretical 
framework to a general problem, that of finding the dis- 
tribution of serving times of individuals forming a crowd 
before a desk. Our model is shown to provide results in 
line with the physical intuition in the case of identical 
agents, which confirms the soundness of our theoretical 
approach. However, we find less obvious results when the 
case of poly-disperse crowds are considered, which con- 
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stitutes a first illustration of the usefulness of our model. 
It should be noted that, while queuing strategy issues 
have been analyzed in some specific contexts of crowd dy- 
namics, such as parallel multi-lane pedestrian traffic [17], 
surprisingly it seems that no simple computational study 
has yet tackled the problem examined in this paper. 

The article is organized as follows. First we describe 
the model and discuss the main results that validate our 
scheme in the case of a homogeneous crowd. Subse- 
quently, we introduce a given degree of poly-dispersity so 
as to examine the statistics of serving times and their de- 
pendence on agents' size. Then, before summarizing and 
discussing our findings, we report an instructive struc- 
tural analysis of the two-dimensional crowd-queuing sys- 
tems. 



II. THE MODEL 

In our model a crowd is represented by a collection of 
hard disks, possibly enclosed in a given perimeter with 
specific boundary conditions. Broadly speaking, the evo- 
lution of the crowd configuration is considered as a two- 
step process. In the first step (i) one agent is chosen 
according to a given rule and is displaced following a 
certain prescription (possibly to infinity, that is removed 
from the game). In the second step (ii) the rest of the 
group is let equilibrate through a Monte Carlo dynamics 
under a pre-assigned displacement rule until a stationary 
state is reached, as measured through the convergence 
of a suitable global order parameter. The details of the 
rules applied in steps (i) and (ii) should reflect the na- 
ture of the process that one wants to simulate. Here we 
introduce a simple set of rules, that could be valid in, or 
else easily modified for, many circumstances where one 
is interested to simulate a flow of people. For example, 
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this could be the case of escape dynamics through doors 
at specific locations, such as in panic stampedes [8, 9] or 
other issues related to the motion of groups of pedestri- 
ans within specified perimeters or across bottlenecks [15]. 

We consider an ensemble of N impenetrable circles in 
the plane of radii r; = 1 + (2z — l)Ar, i = 1, 2, . . . , N, 
where z may be a constant (homogeneous crowd), or else 
be drawn from a set of real numbers specifying the degree 
of heterogeneity. In this paper we either set z = 1/2 or 
consider a uniform distribution in the interval z £ [0, 1]. 
In the specific case chosen to validate and illustrate our 
model, dubbed here crowd- queuing, it should be imag- 
ined that the disks model individuals waiting around a 
counter located at the origin of a Cartesian frame of ref- 
erence and, e.g., selling tickets. The iV-agent crowd is 
served, one individual at a time, following the series of 
serving times t n , n = 1, 2, . . . , N. At each step, the served 
agent is removed from the system, so that the group at 
the n-th step comprises N(n) = N — n individuals. The 
idea is to model a semi-circular geometry, whereby in- 
dividuals disappear as soon as they are served as they 
are admitted beneath the counter (say into the theatre 
or cinema hall), e.g. through a door or down a stair- 
case next to the counter itself. It is clear that, under the 
reasonable hypothesis that walls flanking the counter are 
reflecting (except for the aforementioned exit), it is easier 
to study the equivalent case of a fully circular geometry 
with agents crowding all around the teller and being re- 
moved from the play as soon as they are served [25]. 

For the sake of simplicity, we will take t n+ i — t n = r 
and measure time in units of serving steps. It should 
be noted that this implies a minor loss of generality, as 
we are mainly interested in average quantities, such as 
the average time taken to reach the counter from a given 
distance. More generally, in fact, if P(t) indicates the 
(finite-mean) probability density describing the (uncor- 
rected) intervals between serving times, the average time 
corresponding to m serving steps will still be linear in m, 
namely (At) m = m (t). 

According to our theoretical framework, we approxi- 
mate the dynamics of the group as a sequence of serving 
and rearrangement steps. The basic assumption is that 
the two processes occur on well separated time scales, the 
rearrangement stage being regarded as instantaneous on 
the time-span set by r. At each serving step, the individ- 
ual that is located closer to the counter disappears, and 
the rest of the agents are rearranged. Such process is 
iterated until all agents have been eliminated. Typically, 
we record the serving times of all agents initially belong- 
ing to a certain number of concentric circular shells. By 
doing this, we can measure the average number of serv- 
ing steps as a function of the initial distance from the 
counter, the average being taken over all agents from the 
same shells, as well as over K equivalent initial configu- 
rations of the iV-agent crowd. 

In an ordered, one-dimensional queue the rearrange- 
ment rule is a trivial one - move all agents simultane- 
ously of one position in the direction of the counter. In 



our case, agents have the freedom of moving in two di- 
mensions, whilst they still share a bias towards moving 
along the radial direction, as reaching the counter is ev- 
erybody's ultimate goal. In accordance to such principle, 
we construct the rule for the Monte Carlo rearrangement 
stages as follows. At each MC step j, a radial move is at- 
tempted of the same magnitude Sr(j) for each agent. To 
that, a component 6u(j) is added with probability p along 
a randomly chosen unit vector forming an angle with the 
centripetal direction in the interval [—tt/2, tt/2]. The lat- 
ter moving rule is meant to express the agents' will to 
find alternative routes leading faster to the counter, tak- 
ing advantage of local density fluctuations. The magni- 
tudes of the two types of attempted displacements 5r(j) 
and Su(j) are updated at each step, after the entire set 
of moves is performed, so as to keep the acceptance rate 
a(j) (that is, the fraction of effectively displaced agents) 
around the target value a = 0.5. The acceptance rule for 
each MC move is realized by enforcing hard-core repul- 
sion. 

Following the elimination of the served agent, the MC 
rearrangement is halted when the crowd has reached a 
stable configuration. Obviously, there is no unique man- 
ner to assess whether the system has attained such a sta- 
tus. In fact, also in view of the intrinsic non-equilibrium 
nature of the process that we wish to simulate, such con- 
dition should be understood in a dynamical sense. For 
example, this can be measured by following the evolution 
of a suitable global observable O and fixing the number of 
rearrangement moves M such that \Om — 0a/-i|/Om < 
TOL, where O m = Efeli 0(kAj)/M is the M-step run- 
ning average computed by sampling every Aj MC moves 
and TOL is the required accuracy. In the calibration 
runs, we chose to monitor the convergence of the accep- 
tance rate a and of a global structural indicator, namely 
the overall jamming index T, adapted from Ref. [18] to 
the case of size- heterogeneous disks, 

i=l i=l L ^-yeSsli) J 

where £3(1) is the ensemble comprising the first three 
neighbors of the i-th agent and dij are the center-to- 
center separations. As a results of different calibration 
tests, we found that the choice Aj = 50, TOL = 10~ 4 
yielded a satisfactory compromise between accuracy in 
the convergence of both indicators and simulation speed, 
with M{n)/N(n) ~ const (n being the serving step). For 
the sake of computation celerity, we only monitored the 
observable a for arresting the rearrangement stages in 
production runs. 

The initial configurations for the serving process are 
obtained as follows. The disks are first arranged on 
a square lattice within a circle of radius R = 1 + Ar 
(R = 1 for z = 1/2) and their radii rescaled so as to 
match the required value of the initial disk area frac- 
tion <p = Et=i r i/-^ 2 - A 11 initial cycle of Mi Monte 
Carlo steps is then performed with p = (purely radial 
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p = 0.2 a 



moves) and boundary conditions such that displacements ^ 
are taken modulo 2R. The configurations obtained after 

different scries of Aii ~ O(10 4 ) MC steps are used as p = 0.5 

independent initials conditions for the serving process. 
Of course, arbitrary initial geometries can be enforced in 
the same manner, by retaining only the agents that lie 
within the required perimeter. ^ 

It should be noted that, besides the measure of the K 0.4 

size heterogeneity Ar, our simulation protocol has only 
one parameter, namely the probability p of non-radial 0.2 - 

displacements. 
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In the case of a one-dimensional queue of N individu- 
als, it is clear that an agent initially at a distance d oc n 
from the counter (that is, the n-th agent), will be served 
after a time n (in units of serving steps). In two dimen- 
sions, if the same number of people are crowding around 
the counter, there will be competition among the agents 
lying initially within the shell at distance d from the ori- 
gin. In this case, the equivalent of the one-dimensional 
queuing process will be realized if, from the d— shell, each 
agent will have to wait that all the individuals closer than 
her to the counter be served before she is served. In this 
case, one should expect a number of serving steps given 
by ' 

n seq (d) = N ( ±) (2) 



A. Homogeneous crowds 

The upper plot reported in Fig. 1 shows that, un- 
der the simple rules of our rearrangement dynamics, 
(n(d)) = n scq (d) for a group of identical individuals. This 
means that, as suggested by simple intuition, on average 
there is no difference between an ordered onc-dimcnsional 
queuing process and a chaotic mob formed by an equal 
number of agents crowding around the same serving desk. 
Remarkably, this result appears rather insensitive to vari- 
ations in the probability of non-radial moves. This con- 
firms the soundness and robustness of our approach. In 
the case p = (purely radial moves), we have found that 
the external shells are left virtually motionless, as a re- 
sult of which the crowd is served sequentially from the 
inner shells to the outer ones, more or less following the 
ordered list of distances from the origin in the initial con- 
figuration. Consequently, the sequential serving law (2) 
is (trivially) recovered in the case p = too. Indeed, 
the requirement p ^ reflects a realistic characteristic of 
the dynamics of waiting groups, where individuals share 
an obvious bias towards radial moves that makes them 
approach the counter, but from time to time deviate side- 
ways as local bubbles of low density appear under way. 
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Figure 1: Upper panel: average number of serving steps as 
a function of initial distance from the counter for different 
values of the probability of non-radial moves p. The dashed 
line is a plot of formula (2). Lower panels: distance from the 
counter of some individuals lying initially in a given shell at 
distance do from the origin as a function of the serving steps: 
do = 1 (center) and do = 0.61 (bottom). The arrows indicate 
the average waiting time from the shells. Other parameters 
are: N = 843, p = 0.2, cj> = 0.6, Ar = and K = 30. 



The data reported in the upper panel of Fig. 1 repre- 
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Figure 2: (Color online) Lower left panel: average, maximum and minimum number of serving steps as functions of the initial 
distance from the counter (symbols), n seq (d) (solid line) and n scq (d)/2 (dashed line). Upper panel: (filled squares), 

pq— 0.25 ( enl pty squares), P?. -0 (asterisks) and p^- - 25 (crosses). Right panel: probability density of n(d)/n scq (d) (stairs) and 
normal distribution N([i = 1, a = 0.28) (solid line). Other parameters are: N = 843, <f> = 0.6, p = 0.2, Ar = and K = 30. 



sent averages over many individuals lying initially at the 
same distance from the counter. As a matter of fact, in 
a two-dimensional arrangement different agents from the 
same shell show amply different waiting times, as testi- 
fied by the two cases reported in the two other panels of 
the same figure. It is apparent that serving times up to 
about 30 % below the average characterize certain indi- 
viduals from the same shell, while at the same time others 
will have to wait longer to be served. In order to quan- 
tify the effect of such fluctuations, in Fig. 2 we plot in 
one representative case the average, maximum and min- 
imum serving times as functions of the shell distances 
from the counter. We see that the lowest serving times 
(d) are well approximated by the same sequential 
law (2) but corresponding to half the number of agents. 
Hence, the luckiest individuals arc served as if they were 
crowd-queuing with half the number of people. However, 
the less fortunate agents wait as if they shared the place 
with about twice as much agents, with the obvious con- 
straint that the maximum serving time cannot exceed N. 
It is interesting to remark that the last-of-all condition 
(n max (d) = N) persists for agents initially sitting as close 
to the counter as d = 0.8R. 

Additional, more quantitative information can be ob- 
tained by calculating the probability that an agent be 
served in a number of steps below (or above) a given 
fraction of the sequential-queuing prediction from a given 
distance. Let Afd indicate the ensemble of waiting times 
of all agents from the d-shell over all the independent 
serving runs. Then, we can measure the probability that 
an individual is served in a number of steps lower (or 
greater) than a given fraction 1 ^ q (0 < q < 1) of the 



theoretical expectation, Eq. (2), that is 

P| (d) = «6 [±(1 =F q)n sc(1 (d) T <d)])) n(d)( . Ud (3) 

where &(x) is the Heaviside step function and the opera- 
tion ((...)) makes the double average over shell members 
and independent runs explicit. 

The results reported in Fig. 2 show that the chances 
of reaching the counter in a number of steps below (or 
above) the purely sequential prediction (q = 0) are about 
50 %, and to a large extent insensitive of the initial dis- 
tance. Furthermore, the chances of being served faster 
than 0.75 (or more slowly than 1.25) of such expected 
waiting time are between 15 and 20 %. These observa- 
tions make more quantitative the intuitive feeling that 
competing in a crowd to be served can be rewarding 
as well as imply longer waiting times than in ordered, 
one-dimensional queues. Furthermore, the above results 
suggest that the reduced serving times from all shells, 
n(d) I n seq (d) , are described by the same probability den- 
sity. The right panel in Fig. 2 shows that such curve 
is reasonably well [26] approximated by a normal den- 
sity of unitary mean and standard deviation a = 0.28. 
This provides a handy criterion for quantifying to an ar- 
bitrary confidence level the variability displayed by the 
serving times. For example, in about 70% of the cases, an 
agent starting at a distance d will have to wait between 
0.7n scq (d) and 1.3n scq (d). 

It should be remarked that the above results concern- 
ing minimum and maximum waiting times provide ready 
information under the hypothesis of equally spaced serv- 
ing times t n+ i — t n = t. More generally, the same pre- 
dictions are expected to hold true in the case of uncor- 
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Figure 3: Average number of serving steps as a function of Figure 4: (Color online). Serving time as a function of the 



initial distance from the counter for different values of the 
size heterogeneity amplitude Ar. The dashed line is a plot 
of formula (2). Parameters are: N = 843, p = 0.2, 4> — 0.6, 
K = 30. 



agent size for individuals lying initially at a distance d > 0.9R. 
The serving times have been binned along the radius axis and 
normalized to the average serving time of the considered shell. 
Solid lines are guides to the eye. Parameters are N — 475, 
p = 0.3, Ar = 0.3 and K = 30. 



related, normally distributed serving times (likely a rea- 
sonable description of most common counter dynamics), 
provided the serving and rearrangement time scales are 
still uncoupled. In situations where different instances of 
the distribution V(t) are relevant, the relation between 
serving times and number of serving steps could be non- 
linear. In such cases, the convolution of the queuing and 
the serving time statistics should be accounted for explic- 
itly and our algorithm should be modified accordingly. 



B. Heterogeneous crowds 

We have seen that, on average, an individual involved 
in a queuing process within an homogeneous crowd gets 
to the counter following the sequential rule, that is after 
all agents closer than her to the origin are served. Fig. 3 
proves that such average behavior persists also in the case 
of a poly-disperse group of agents, no matter the spread 
of the size distribution. However, each point in Fig. 3 is 
worked out by averaging over a population of agents that 
share the same initial distance from the counter but are 
heterogenous in size. What happens if all those serving 
times are sorted as functions of the agents' radii? 

Fig. 4 provides a clear-cut answer in the case of indi- 
viduals initially lying within the outermost shell: start- 
ing from the same distance, the smallest agents get to 
the counter about 10 % faster than the bulkiest ones. 
This phenomenon seems to depend to a certain extent 
on the details of the rearrangement dynamics, although 
the overall advantage displayed by the tiniest agents over 
larger ones seems not to be altered by the actual value of 
the probability of non-radial moves p. We remark that, 
as expected, the average serving time that characterizes 
the shell (indicated as (n) in the figure) also marks the 



dynamics of individuals of average size, that is those with 
r = 1. 



IV. STRUCTURAL PROPERTIES OF THE 
HARD-DISK CROWD 

It is interesting to analyze quantitatively the spa- 
tial configurations describing our granular crowd systems 
during their temporal evolution. Fig. 5 shows four snap- 
shots of the group at different stages for an initial disk 
area fraction (f> = 0.6. Overall, a tendency of the group 
of agents to deviate from the circular shape is observed, 
likely as a result of non-radial moves performed along 
the way. A similar scenario is confirmed by many other 
analyzed instances (not reported here). Likewise, it ap- 
pears that the serving process also leads to a tightening 
of the system. In order to quantify this effect, we have 
calculated the variations of the agents' effective surface 
coverage by performing Dirichlct tessclation of the con- 
figurations at the end of each rearrangement stage. Fig. 6 
reveals that the first phase of the process indeed marks a 
rapid increase of <j> between 10 and 20 %, depending on 
the initial value of the disk area fraction. Subsequently, 
the serving process leads to a loosening of the group, 
which appears to be more marked the denser the starting 
configuration. Interestingly, towards the end of the pro- 
cess the last-served agents seem to organize themselves 
into tighter structures. 

Visual inspection of the systems' patterns suggest that 
local order appears to be enhanced late in the serving pro- 
cess (see again Fig. 5). To make this observation more 
quantitative, we have calculated the radial distribution 
function g2(r) (RDF) at different steps in the course of 



Figure 5: (Color online). Four snapshots during the crowd- 
queuing process at different steps n. From top to bottom 
and left to right: initial configuration, n = 300, n = 500 and 
n = 700. The long- and short-dashed lines mark two shells 
at d — 0.99 and d — 0.6, while the corresponding agents 
are represented by dark and light blue circles, respectively. 
Parameters are: N = 843, (f> = 0.6, p = 0.2 and Ar = 0. 
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Figure 6: Effective disk area fraction as a function of the 
serving step during the crowd-queuing process. Parameters 
are p = 0.2, Ar = and K = 30. 



the serving dynamics. The results are summarized in 
Fig. 7 for one typical case. The initial configurations, 
as one should expect, reflect the known RDF of hard 
disk fluids [19, 20]. However, as the serving process pro- 
gresses, short-range translational order gets manifestly 
enhanced alongside with the rapid increase of surface cov- 
erage. 

It is interesting to observe that the calculated RDFs 



Figure 7: (Color online) Radial distribution function versus 
pair- wise distance in units of the diameter D at different steps 
of the serving process Other parameters are: N = 843, (f> = 
0.6, p = 0.2, Ar = and K = 30. 



appear remarkably similar to the ones characterizing 
the non-equilibrium configurations obtained in Ref. [21] 
through a modification of Eden's algorithm targeting a 
pre-determined degree of oricntational order. In fact, it is 
instructive to investigate whether the observed enhance- 
ment of translational order is also accompanied by bond- 
orientational order. To do this, we employed the com- 
puted Dirichlet tesselations to evaluate the global bond- 
orientational order parameter 
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Figure 8: Global bond-orientational order parameter ver- 
sus fraction of removed geometrical bonds during the crowd- 
queuing process (proceeding left to right) for different choices 
of N and <f>. Other parameters are p — 0.2, Ar = and 
K = 30. 
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1 n* 

^ = ^EE cos ( 6 %) ( 4 ) 

b i 3=1 

where rij is the number of geometrical neighbors of the i- 
th disk, 9ij is the angle formed by the bond between the i- 
th and j-th disks in the dual Delaunay triangulation and 
some arbitrary but fixed direction (in our case the x axis) 
and Nb is the total number of such geometrical bonds. By 
definition, ^6 — 1 m the most ordered arrangement, that 

— 1/2 

is the triangular lattice, whereas ij) oc N b in a dis- 
ordered, gas-like phase [27]. Fig. 8 shows that queuing 
processes that start at moderate surface coverage do not 

— 1/2 

develop appreciable bond-ordering, as ip oc N b for the 
whole duration of the process. When the initial area frac- 
tion is higher, the first configurations indeed feature non- 
negligible bond-orientational order. However, the serving 
dynamics invariably brings the system back to gas-phase 
configurations on a rather short time scale - when about 
20 % of the agents have been served, any traces of bond- 
orientational order have been already wiped out. 

Overall, the above analysis reveals that the configu- 
rations obtained through our simple method are typical 
non-equilibrium configurations of hard disks, featuring 
rather high densities and maximum degree of disorder 
at the same time. Therefore, our method might as well 
prove interesting for non-equilibrium hard-disk packing 
issues [21]. 

V. CONCLUSIONS 

In this paper we have introduced a simple agent-based 
Monte Carlo model of crowd dynamics, with the aim of 
putting forward a viable and handy alternative to more 
detailed and involved modeling strategies for a variety 
of crowd dynamics issues. Agents are modeled as hard 
disks and the evolution of the crowd spatial configura- 
tion is reconstructed through a series of two-step moves: 
in the first step a single agent is selected and displaced 
according to a certain rule, possibly to infinity, which 
amounts to removing her from the game (Rl). In the 
second step, the rest of the group is let rearrange through 
a sequence of Monte Carlo moves obeying a certain dis- 
placement rule (R2) until the convergence of a suitable 
global observable has been achieved. The two aforemen- 
tioned rules Rl and R2 are meant to identify the nature 
of the specific crowd dynamics that one wishes to simu- 
late. 

We have shown that out theoretical framework pro- 
vides results in agreement with physical intuition in the 
problem of crowd- queuing, namely that of a group of 
agents who gather in a disordered manner around a 
counter (a ticket-desk, a bar, etc.) waiting to be served 
one by one and leaving the place once served. Our 
model shows that on average the crowd reproduces one- 



dimensional, sequential queuing for a group of agents of 
identical size, in agreement with intuition. This means 
that an agent starting from a certain distance d from 
the counter will have to wait that all agents initially en- 
closed within the circular d-shell will have been served. 
Obviously, however, a two-dimensional crowd provides a 
queuing environment that can yield slower-than as well 
as faster-than-sequcntial serving times for a given indi- 
vidual. Our results indicate that the two instances occur 
with a probability of about 1/2, independently of the ini- 
tial agents' position. Remarkably, the above results are 
rather insensitive to changes in the probability of non- 
radial moves, which testifies to the robustness of our ap- 
proach. 

Remarkably, size of agents matters in crowd dynamics 
processes. On average the serving times of an hetero- 
geneous group still follow the intuitive sequential predic- 
tion, independently of the strength of the agents' size het- 
erogeneity. However, starting from identical initial dis- 
tances, the smaller the agent the shorter her waiting time. 
The tiniest individuals manage to sneak through their fel- 
low queuers more effectively, thus reaching the counter 
the first. It is intriguing to draw an analogy between 
our findings, i.e. larger disks lagging behind smaller 
ones, and the well-known Brazil-nut effect, whereby the 
largest particles migrate towards the surface when a gran- 
ular heterogeneous mixture is shaken [22]. In fact, par- 
ticle size segregation effects have already been reported 
in two-dimensional hard-disk packings [23]. In the case 
of shaken mixtures, the bias provided by gravity plays 
an important role in determining the segregation phe- 
nomenon. Within the framework of our rearrangement 
dynamics, the agents' mixture is Monte-Carlo shaken un- 
der the bias of preferred origin-pointing displacements. 
In this sense, explanations of the Brazil-nut effect that 
focus on infiltration of small particles into voids created 
underneath larger ones during shaking might provide rel- 
evant clues to our process. Conversely, the analogy with 
crowd dynamics might prove useful in the rationaliza- 
tion of size effects in poly-disperse mixtures. For exam- 
ple, this could be the case of the emergence of indirect 
attractive forces favoring flocculation of larger particles 
resulting from the interference between topological per- 
turbations (depletion wakes) induced by the large parti- 
cles [24]. 

Finally, we have performed a structural analysis of 
crowding ensembles. We have shown that the rearrange- 
ment dynamics results in an overall increase of the agent 
area fraction, enhancing at the same time short-range 
translational order in the system as the serving process 
progresses. On the contrary, orientational order is effec- 
tively suppressed, even in the case of high initial packing, 
leading to rather dense and at the same gas-like, highly 
disordered non-equilibrium configurations, reminiscent of 
those accessible through certain seed-based growth pro- 
tocols [21]. 
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